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Abstract 

In this paper, a probabilistic representation of the tau functions 
of KP (Kadomtsev-PetviashviH) sohtons in terms of stochastic areas 
will be presented. 

1 Introduction 

In the introduction, after giving a very short introduction to the theory 
of solitons following [16], we recall some existing results from probabilistic 
approaches. 

1.1 Solitons, tau- functions, and Sato's Grassmannian 

By solitons, we usually mean solitary wave solutions (behaving like a particle) 
to a class of non-linear wave equations including the KdV (Korteweg-de Vries) 
equation 

du 1 d^u 3 du , . 

'dt ^ ia^^ 2^9^ ^ ■ 

as the most notable example. 

The first giant step in the study of solitons was made by Gardner, Greene, 
Kruskal and Miura [3], where they observed that (i) the eigenvalues of the 
Shorodinger operator 
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where m is a solution to f ll.ip . is constant in time parameter t, and (ii) one 
can construct a soliton solution to f ll.ip by applying the inverse scattering 
method, by which we mean the (mathematical) method to construct (un- 
known) potentials out of given scattering data, which had already been fully 
developed. The relation is most clearly seen when the potential is reflection- 
less as 

u{t,x) = 2— \ogdet{I + G{x,t)), (1.2) 



where 



G{x,t) 



mimjC 



Vi + Vj 

The constants rij,mj, j = 1, ■ ■ ■ ,n are so-called scattering data. 

The observation (i) together with the awareness of the existence of the 
infinite invariants in [3] motivated another seminal paper by P. Lax [13] . 
where the KdV equation f ll.ip is understood as the compatibility between 
the two equations: 

+ u{t, x) ) w{=: Pw) = Kw, [k is an eigenvalue) 



This compatibility is rephrased as the celebrated "Lax equation" : 

^ + [P,B] = 0, (1.3) 

where the bracket is the commutator; [P, B] = PB — BP. 

By considering pseudo differential operators such as for G N and 
their infinite series, we have in fact B = (P^/^)_|_, where (-D)+ is the differ- 
ential operator part of the pseudo differential operator D. In this Lax form, 
the existence of the infinite many invariants can be rephrased as 

dP 

— + [P, (P'^/^)^] = 1, 3, 5, ■ ■ ■ , 2n + 1, ■ • • , 

where u = u{xi, x^, - ■ ■ , X2n+i, ■ ■ a function of infinitely many variables. 
The KdV case fll.3p is retrieved by setting xi = t,x^ = x. Each Lax equa- 
tion generates a non-linear evolution equation with respect to X2k+i since 
[P, (P'^/^)_|_]'s are all multiplication operators. The totality of the generated 
equations is usually called KdV hierarchy. 
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If we instead start with the operator 



L = d + J2^j9-^^ 
j=i 

then we still have that [L, {L'')^] are all multiplication operators, and hence 
we obtain infinitely many nonlinear differential equations with respect to m/s 
of infinitely many variables xi,X2, - ■ ■ ,Xn,--- by the Lax equations: 

— + [mL%]=0, A; = l,2,---. 

The family is called KP hierarchy since the KP (Kadomtsev-Petviashvili) 
equation, 

Sd'^ui d f dui 3 dui Id^ui 
4 dxl dxi \dx3 2^^dxi 4 dxl 

which is easily seen to be a generalization of the KdV to a two dimensional 
model, is deduced from the equations with k = 2 and k = 3. The KP 
hierarchy as a whole is also a generalization of the KdV hierarchy since the 
latter hierarchy is obtained by a reduction {L"^)- = from the former. 

The equations in KP/KdV hierarchy are all "soliton equations" in the 
sense that they all have exact solutions of soliton type0. In fact, according to 
Sato's theory of infinite dimensional Grassmannian ([I9], see also [I6l[18]), all 
the Uj 's of the hierarchy are simultaneously generated from a single function 
called tau-function t in the following way: determine wi, W2, etc, by 

rjXi- l,X2- ^,---) _^ ^ Wi ^ W2 ^ 

t{xi,X2,---) k k'^ 

by comparing the coefficients of k~^ , j = 1, 2, ■ ■ ■ , and then ui, U2, etc by 



For example, we have 



' dx\ 



«i = 27z:ilogr. (1.6) 



-•^The solitons are basically rational functions of the exponential functions of the form 
^cijXj j-Qj. gQjj^g constants c^j 's. 
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In particular, we see that if r is a polynomial of e^'^'J^-''s, then u/s are all 
"solitons" in that they are all rational functions of e^^'^^^ 's. 

The tau functions are characterized as a solution to a family of quadratic 
differential equations called Hirota equations, which are nothing but Pliicker 
relations that define Sato's infinite dimensional Grassmannian. That is to 
say, a tau function of the KP hierarchy is a point in the Sato's Grassmannian. 
It should be noted that in the Sato's theory, the KP hierarchy is the most 
universal one, out of which many well-known soliton equations are obtained 
by a reduction. 

The following functions are known to be among the tau functions of the 
soliton solution of the KP equation: 

r(xi,X2, ■ ■ ■) 

for / = {1, ■ ■ ■ ,n}, n G N, where mi, ■ ■ ■ ,m„, pi, ■ ■ ■ ,p„, and Qi, ■ ■ ■ ,qn are 
(indefinite) constants. This is alternatively written as 

t{xi,X2, ■ ■ ■) = det(/ + G{xi, X2, ■ ■ ■ ))> (1-7) 

where 

\Pi-(lj Ji,j 

The formula f ll.7p is a generalization of f ll.2p since we retrieve it by f ll.6p and 
the reductions of qj = —pj, = for Z > 4. 

Remark 1. It should be noted that, if / is a solution to a Hirota equation then 
so is Ce^J '^j^j for arbitrary constants C, ci, C2, ■ ■ ■ . Therefore tau function 
is stable under the multiplication of the factor Ce^^'^^^^. This property will 
be used in the proof of Theorem [31 

1.2 Probabilistic approach to solitons 

As far as we know, the first attempt to represent solitons in terms of the 
expectation of Wiener functionals was made by S. Kotani [11] in 2000. Ac- 
cording to [8], Kotani constructed the following correspondences. Let S be 
the set of all pairs ((T+,cr_) = a of non- negative measures each on R_ such 
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that Jj^ 



-lA 



<y±{dX) < oo for any A > 0. For a G S, associate a Gaussian 



process with mean whose covariance C{u,v 
by 

1 



C(m, v; a) 



4 

1 

+ 4 



-1/2 



/— z(m+i') 



— z|m— i"! 



E[X{u)X{v)] is given 
cr+((iz) 



-1/2 



R 



-z(u+v) 



Let Q be the totahty of the function q'^ with a G S, where 

1 



g [X) 



rj2 

-4— logE 



exp 



\X-iy)\'ds 



:i.8) 



Then Kotani showed that Q is the closure (with respect to the topology 
of uniform convergence on compacts) of UQ„, where Q„ be totality of the 
reflect ionless potentials of scattering data consisting of 2n constants. 

In a somewhat different line, K. Hara and N. Ikeda [5] derived from the 
Fourier transform of a class of quadratic Wiener functionals a dynamics in 
the Grassmannian as a finite dimensional analogue to the Sato's framework 
dH-dlS]) etc. 

Soon after that N. Ikeda and S. Taniguchi [8J obtained a specific and 
more "stochastic analysis oriented" construction of the Gaussian process X'^ 
in (11. 8p than Kotani's method. They set 



:i.9) 



where a > 0, c G R" , p G R" and is an Ornstein-Uhlenbeck process in R"' 
starting at defined as the soliton to the following SDE: 

d^t = dWt + diag{pi, ■ ■ ■ ,Pn}^tdt. (1.10) 

The measure a in Kotani's correspondence is given as 

(7±{du) = cf6_p2{du). 

Related studies and surveys concerning the quadratic Wiener functionals can 
be found in [20H23] . and more recently in [9]. 

Here we remark that all the probabilistic results cited here are on KdV 
solitons, and not extendable to KP. In this paper, we will present a proba- 
bilistic representation of KP solitons using generalized stochastic areas (see 
Theorem [2] and Theorem [3]) . 
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1.3 Organization of the present paper 

In section [21 we will introduce Levy's stochastic area formula and present its 
generalization as Theorem [2] and its proof. Then in section [3l we will show 
that the generalized stochastic area formula is parameterized as a tau func- 
tion of KP solitons. In section HJ we will give a probabilistic interpretation of 
the reduction from KP- to KdV-solitons. In section we will present another 
realization, where the dimension of the Wiener space is fixed to two. 



2 A generalization of Levy's stochastic area 
formula 

Let {fl, J-", P) be a probability space and W = (W^, W^) be a two-dimensional 
Brownian motion on it. The area enclosed by the curve s t— )■ Ws and its chord 
up to time t, which is usually called stochastic area of W, is given by 

St ■■= \ dwl - £ wl dwi\^ . 

The characteristic function of St is explicitly given as 

E[ev^«^*] = (cosh|)-i (eeR), (2.1) 

and conditioned one is also given explicitly as 

E[e^^^*|iy/=x,iy,^ = l/]= cotha) (^^j^)^ (2.2) 

2 sinh ^ 

which were found by Paul Levy [i4j using Fourier series expansion of W . 
Either is often called Levy's (stochastic area) formula(s). There have been 
plenty of studies related to the formulas. For example, the heat kernel of 
the Heisenberg group can be obtained by a slight modification of the formula 
([1], see also (ID])- Many alternative proofs and generalizations have been 
found (121], P, 0, [6], [7], etc). 

In this paper, we give the following generalization of (12. ip . In its proof, 
the second Levy formula (12. 2 p plays a crucial role. 
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Theorem 2. Let = (VT''^, VT''^), I = 1, ■ ■ ■ ,n be mutually independent 
two-dimensional Brownian motions starting at the origin, and stochastic ar- 
eas of will he denoted by 



Jo 



Let A := diag{Ai, ■ ■ ■ ,Xn}, where Xi,l = 1,2,- ■■ ,n are positive numbers. 
Let A = {o,ij)i<ij<n be a real n x n matrix, and be its symmetric and 
skew- symmetric part respectively, namely, = {A± A*)/2. Denote = 
{Wl'\ ■ ■ ■ , ly"'^) for 1 = 1,2, and define for z e C 

S{z) = SaA^) 

n 2 
1=1 i=l,2 

(2.3) 

Then, if either maxilXil or \\C^\\ is sufficiently small, we have 

cosh Ai + ai^i sinh Ai ai^2 sinh A2 ■ ■ ■ ai,n sinh A^ ' 

02,1 sinh Ai cosh A2 + 02,2 sinh A2 ■ ■ ■ a2,n sinh A„ 

a„_i sinh Ai a„^2 sinh A2 • • ■ cosh A^ + a„^„ sinh A^ 

Proof. We first calculate the conditional expectation of e'^*^""-' conditioned by 
Wi = (W{, Wj) . By the Levy's formula (12. 2p with some analytic contin- 
uation, we have for sufficiently small o" G R (such that the random variable 
gO-j]A;S jg integrable), 

= -. r^^P ^ {aXi cot aXi - 1) . 

sm aXi V 2 / 



Therefore we have 

(tXi_ ( 1 



E[e^W|Wi] = TT exp ( - ^((M(a) - I + C(a))Wi, Wi)) 

sm cA; \ A J 



I 
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where 

^ /crAcotcrA 

^^^^=[ aAcotaA 
with 

cot crA := diag{cot aAi, ■ ■ ■ , cot aXn} 



as usual, and 



-aA^C-Al (T^A^C+A^ 



Since ||M(cr) + C{a) — I\\ — t- as cr — )■ 0, we can take a small enough 
to ensure that M(cr) + C(cr) is positive definite. Then, applying quadratic 
Gaussian formula for such a, we obtain 

E[e^('^)] = n det(M(cr) + C(a))-i (2.5) 

sin cr A; 

1=1 ' 

We may go further as 
det(M(a) + C((t)) 

V -(tA^C^As (tA5(A5 cotaA + (tC+A^) ) 

= det((TA5(cot aA + aC^ + v^C")A5) det((TA5(cot aA + aC+ - v^C")A^ 
(Since C~ is skew symmetric) 

YliaXi) det(cot aA + (tC+ + v^C" 

Hence ( 12. 5 p is turned into the following equahty: 

E[e^(")] = det cos aA + {(tC+ + v^C") sin ctA) ^ (2.6) 



where sinaA := diagjsin cxAi, ■ ■ ■ , sincrA^}. 

The right-hand-side of f l2.6p is apparently meromorphic in a G C. Now, 
we want to see if an analytic continuation to a domain including z = 
is possible or not. To see this, it suffices to check the differentiability of 
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the left-hand-side of fl2.6p with respect to a. Namely, we need to check the 
integrability of 

dz 



E 



. i=l 1=1,2 



Since 5* is quadratic Gaussian, the integrability is inherited from that of e'^'-^-' 
itself, which is guaranteed if either max/ \ \ or ||C^|| is sufficiently small. □ 



3 Parameterization as a tau function of KP 
solitons 

As we have stated, a tau function r of the n-soliton solution of the Kadomtsev- 
Petviashvili equation (KP equation) is expressed by 

r(xi,X2, ■ ■ ■ ) = det(/ + G'(xi,X2, ■■■))) (3-1) 

with 

V Vi Qj J l<ij<n 

where 

6 = {Pi - + {Pi - qf)x2 + ■ ■ ■ , i = 1, ■ ■ ■ , n. 



and rrii > 0, pi and qt are parameters. 

Theorem 3. Let P = ( — ^— ) , and assume that miuj ,,- \pi — qj\ is 

\Pi 1j / l<ij<n 

sufficiently large so that I + P is invertible. Then, if we put A = {I — P){I + 
P)~^ and A := diag{ — + logmi), ■ ■ ■ , — |(^n + logm„)}, we have that 
(E[e'^'^'*''^^''])~"^, where Sa,a is defined by l{2. 3\) . defines a tau function of 
KP solitons. 

Proof. Since 

G = e-^Pe-^, 
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we have 



r = det(/ + e-^Pe-^) 
= det e"^ det(e^ + Pe'^) = det(/ + Pe^^^) 



On the other hand, 



det (cosh A + A sinh A) 



det + A- 



2 2 

= 2-" det { (/ + A)e^ + (/ - A)e-^} 

= 2-" det{(/ + A)e^} det (/ + (/ + - A)e-^^] 

= 2-" det(J + A)e-^ E(e.+iogm.) j^j ^ pg"^^) . 

The last equahty follows since 

A = {I - P){I + P)-^ ^ P = {I + A)-\I - A). 

As we have stated in Remark^ 2~" det(/ + A)e~5(^'+'°§'"') is a trivial factor 
and thus by Theorem [2] we have the assertion. □ 

4 Reduction to Ikeda-Taniguchi's construc- 
tion 

As we have discussed in section [TTT| we have f ll.2p by the reduction of qj = —pj 
in (11. 7p . In this section, we review this from the perspective of stochastic 
analysis. We will show that when C~ = 0, the expectation of the exponential 
of the generalized stochastic area is reduced to that of the exponential of 
the time integral of an Ornstein-Uhlenbeck process, which corresponds to 
the Taniguchi-Ikeda's construction (II. 8p . (II. 9p and (ll.lOp of reflectionless 
potentials/tau functions of KdV solitons. 
Precisely speaking, we have the following 

Proposition 4. Suppose that A in Theorem\^is symmetric. Then 
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where X^'^ = ((A — AKA)^^^) and ^ is an Ornstein-Uhlenbeck process on 
H'^ starting at and satisfying 

d^t = A^dBt + AA^t dt, (4.1) 

with B being an n- dimensional standard Brownian motion. 

Proof. We first note the following identity since its right-hand-side also equals 
to that of fl2.4p with a replaced by (see e.g. pSj): 

Then since = A and C~ = 0, we have 

1=1,2 

'E[e-^'i^o(w's'')^dsmAiAAhwl,wl)^y _ 

By applying Ito's formula, 

e" 4- loiwt'? ds+liAiAA^wlwl) 

^ ^^trAA^J^iA^ AAiwld'^D-y^ \Ai AA^Wll^ ds^-y^{{A-AAA)Ahwl,Ahwl) ds _ 

Define Q by 



dQ 
dP 



= e 



^ll,AhAA\-N\,d-N\)~\\l\AhAA\^\\^ds 



Then by the Maruyama-Girsanov theorem, we see that W under Q has the 
same law as ^ of (14.11) . This completes the proof. □ 



5 A realization in 2D- Wiener space 

In this section, we shall show that the 2r;,-dimensional Brownian motion used 
to represent ra-solitons in Theorem |2] can be replaced by a 2-dimensional one 
irrespective of n. This observation suggests that the whole Sato's Grassman- 
nian would be realized in 2-dimensional Wiener space. 
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Let W = {W^, W'^) be a 2-dimensional Brownian motion starting at the 
origin, {e^ — {e\, • ■ ■ el^) , i — 1, ■ • ■ n} he an orthonormal basis of R", and set 

n 

fi{t) := \/n Ve;lri^ i = 1,2,--- ,n. 



i=\ 



Define 



and 



:= ^ fAs) dWl^ Mt) dWl -I Ms) dW^^ m dWl 

In this section we assume that Aj > for all i. We shall denote the 
entry of the matrices AzC'^'Aa and A2(/ + C~)A2 by A^^, and A^^-, 
respectively. Note that A^ — \. We assume that either max; |A;| or ||C+|| is 
sufficiently small to ensure the integrability. 

Proposition 5. We have that 

= det(coshA + AsinhA)"^ 
Proof. We will show the following equivalence in law: 



\i=l i=l,2 / 

Let us start with the following direct calculation: 

a=l,2 k,l 

and 

k<l 
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By the scaling property of the Brownian motion, the process 

{(W^(U/-l)/n - ^ll-D/n. K^l-l)/n ' ^(i-l)/n) : < 5 < 1, / = 1, • ■ ■ , n} 

is identically distributed as 

{n-^'^{Wl^\ : < s < 1, / = 1, ■ ■ ■ , n} 

Here {VF''^, M^''^, I — 1, - ■ ■ ,n\ are 2n-dimensional Brownian motions staring 
at the origin. In particular, 

I r {W-Wl_,y^)dWl-{Wl-Wli_,y^)dWl}^ = lr-- ,n\ 

yj{l-l)/n J 

is identically distributed as 



n 



J^' (wl^'dWt' - Wt'dWl'') , / = 1, • • • , n| , 



and 

{(W^/n - Wl,_,yJ{W^/^ - W^(ti)/J, l<k,l<n} 

is identically distributed as {V^f '^VF^I'^, 1 < k,l < n}. Therefore, we have the 
following identity in law: 

i,j a=l,2 k,l i,j 



1=1,2 



a=l,2 
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and 



k<l i,j i,j 



= fiT^ATWldW])- [\T*ATWldW^) (5.2) 
Jo Jo 

+ (T*Aic-A3TW},W?) 

n „i 

4 5^ A J (wt'dWt'-Wt'dWt^) 
1=1 



Here W", a = 1, 2 are n- dimensional Brownian motions, and T := [e^ ■ ■ ■ e"], 
which is an orthogonal matrix since {e^ ■ ■ ■ e"} is an orthonormal basis. Note 
that A2C~A2 is also skew-symmetric and thus (AaC^Aa)^^ = for i = 
1, ■ ■ ■ ,n. 

Using the equivalences (15.11) and (15. 2p in law, we can establish 

E[gE,,i(v^A^,j5,,,+|A+.5+ )j ^ ^ det(cosh A + A^AA^ sinh A)"^ 

in the same manner as we have done in the proof of Theorem O Since A 5 
and sinhA commute, we have the assertion. □ 
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